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We investigate the effect of friction on domain formation in disoriented chiral condensate. We 
solve the equation of motion of the linear sigma model, in the Hartree approximation, including a 
friction and a white noise term. For quenched initial condition, we find that even in presence of 
noise and dissipation domain like structure emerges after a few fermi of evolution. Domain size as 
large as 5 fm can be formed. 
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The possibility of forming disoriented chiral condensate (DCC) in relativistic heavy ion collisions has generated 
considerable research activities in recent years. The idea was first proposed by Rajagopal and Wilczek They 
argued that for a second order chiral phase transition, the chiral condensate can become temporarily disoriented in 
the nonequilibrium conditions encountered in heavy ion collisions. As the temperature drops below T c , the chiral 
symmetry begins to break by developing domains in which the chiral field is misaligned from its true vacuum value. 
The misaligned condensate has the same quark content and quantum numbers as do pions and essentially constitute 
a classical pion field. The system will finally relaxes to the true vacuum and in the process can emit coherent pions. 
[~^. ' Since the disoriented domains have well defined isospin orientation, the associated pions can exhibit novel centaurodike 
t— I , |^| H fluctuations of neutral and charged pions [p|-fi2| . 

Most dynamical studies of DCC have been based on the linear sigma model, in which the chiral degrees of freedom 
are described by the real 0(4) field $ = (a, II), having the equation of motion, 
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The parameters of the model can be fixed by specifying the pion decay constant, f n =92 MeV and the meson masses, 
m,r=135 MeV and m (T =600 MeV, leading to A = (m 2 -m 2 )/2/ 2 =20.14 and v = [(ml-3ml)/(ml-ml)]^ 2 f n =86.71 
ON ' MeV and H = (120.55MeU) 3 fL3fl . It is apparent from eq.[j]that the vacuum is aligned in the a direction & vac — (fn, 0) 
and at low temperature the fluctuations represent nearly free a and 7r mesons. At very high temperature well above 
v, the field fluctuations are centered near zero and approximate 0(4) symmetry prevails. 
Q_i' It is instructive to decompose the chiral field, 



$(r, t) =< <j>{r, t) > +5(f>{r, t) (2) 



where < >is the mean field and S(j) are the semiclassical fluctuations around < <f> > and can be identified with 
quasi-particle excitations. Using eq.^| and taking the average of eq.[l], the equation of motion for the mean fields in 
\ the Hartree approximation can be obtained as p3| , pd[ , 

9 ^/ > - V 2 < <t> >= A(w 2 - <(/) > 2 -3 < d(/)j >-< S<Pl >)<(/) > +Hn a (3) 

where < <f> > 2 ~< cf>i >< <f>i >, S(f>\\ is the component of the fluctuation parallel to < cf> > and S<fi± is the orthogonal 
component. This equation imply that the motion of the mean field is determined by the effective potential, 

V{< (f>>) = j(< (p > 2 +3 < S(f>l > + < S<Pl > -v 2 ) (4) 

which clearly differs from the zero temperature one in presence of fluctuations. By varying the fluctuations, chiral 
symmetry can be restored or spontaneously broken. It is also evident that the evolution of the mean field critically 
depends on the initial values of the fluctuations. When S 2 = (3 < Sip 2 > — < 84)\ >)/6 is large enough the chiral 
symmetry is approximately (as 0) restored. If the explicit chiral symmetry breaking term is neglected, the phase 
transition takes place at the critical fluctuations S 2 = v 2 /6. For 5 2 < 8 2 , the effective potential takes its minimum 
value at < 4>>— (tr e ,0), where a e depends on S 2 . When the mean fields are displaced from this equilibrium point to 
the central lump of the Mexican hat (< <fi >^ 0), the effective mass square 
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will become negative and DCC can form. Since the domain size is directly related to the time scale, during which the 
effective mass remains negative, it strongly depends on the initial condition of the system. By varying the < <pi > 
and S 2 , quench or annealing like initial condition can be obtained Eifl , In an important paper Asakawa et al El] 
studied eq.0 with initial condition corresponding to quench and annealing. They found that domains of disoriented 
chiral condensate with 4-5 fm in size can from through a quench. Annealing on the otherhand, leads to smaller sized 
domains. 

In the present paper, our interest is to investigate the effect of friction on DCC domain formation. Dissipative 
effects like friction damp the motion of the fields, inhibiting large oscillations. Moreover, fluctuations-dissipations 
theorem require that dissipation be associated with noise. One would then expect large reduction in the DCC domain 
formation if friction is present in the system. This naive expectation was found to be true by Biro and Greiner 
p5[ . Using the Langevin equation for the linear sigma model, they have investigated the interplay of friction and 
white noise on the evolution of the order parameter. While noise greatly diminishes the possibility of DCC domain 
formation, in some orbits, large instabilities can result, producing DCC domains. We have also studied the effect 
of friction on DCC domain formation using the Langevin equation JTr| . There we found that for one-dimensional 
expansion on average large DCC domain can not be formed. However, in some particular orbit large instabilities can 
occur. This possibility also reduces with introduction of friction. However, if the friction is large, the system may be 
overdamped and then there is a possibility of DCC domain like formation. Present paper is an extension of the above 
study, the spatial part, which had been integrated out in our earlier analysis is being studied here. 

Appropriate coordinates for heavy ion collisions are the proper time (r) and the rapidity (Y). The change can be 
effected by the following replacement , 

dp ~ Ih 2 7a7 r 97 ~ T 2 W 2 ^' 

It can be seen from the above equation that with the introduction of proper time and rapidity, a dissipative term 
comes into effect in the equation of motion. To illustrate the role of friction, we further introduce a dissipative term (77) 
in the equation. However, fluctuations-dissipation theorem require that dissipation be associated with fluctuations. 
We thus include a white noise term also. To simplify our calculation, we assume boost-invariance in the system, Eq.|| 
can be written as, 

d 2 < > A s d < <f> > d 2 < <j) > d 2 < <\> > w 2 , , „ ln . , 
<^f h + ^ 1Fr~ = &^ ' dy^ h ^ <4>> -T-/2] <(j>> 

+Hn a + ((T,x,y) (7) 

where rj is the friction coefficient. £(t, x, y) is a Gaussian noise of temperature T with correlations, 

< Q{T,x,y) >=0 (8a) 
< Co(ti, xi, yi)(b{T2, x 2 , y 2 ) >= 2i]TS(ti - t 2 )S(xi - x 2 )S(yi - y 2 )S a b (8b) 

It may also be noted that we have replaced the fluctuation term 3 < 6(f)?, > — < 8if>\ > by its counterpart (T 2 /2) 
in the finite temperature field theory 0,|l8| . The equation of motion of fields then depend sensitively on the initial 
temperature of the system. In the following, we assume initial temperature to be T c — \/2f 2 — 2m 2 /X=123 MeV at 
the initial time 73= 1 fm. The cooling of the system is described by the following equation appropriate for scaling 
expansion in 2-dimension, 

T 2 

The friction coefficient was assumed to be rj = rj^ + rj a , where r) Kja - are the friction coefficient of the pion and the 
sigma fields. They have been calculated by Rischke pjj , 

,4A/^ 2 ml r Am 2 1 - expj-m^/T) 1 

l _ ! ' m 2 1 — exp(— m 2 /2mT T T) exp(m 2 — 2m 2 )/2m p iT) — 1 
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We solve the set of partial differential equations |?] with the quenched initial condition. Accordingly the initial fields 
are randomly distributed to a Gaussian form with the following parameters, 
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< a >= (1 - f(r))U (11a) 
<7r,>=0 (lib) 
< a 2 > - < a > 2 =< ir 2 > - < tt, > 2 = v 2 /4f(r) (11c) 

<(7>=<7f J >=0 (lid) 

< g 2 >=< 7T > 2 = u 2 (lie) 

The interpolation function 

f(r) = [l + exp(r-r )/T)}- 1 (12) 

separates the central region where the initial field configuration is different from their vacuum expectation value. We 
use r =5 fm and T=0.5 fm. 

The equation of motion was solved for 500 trajectories. For each trajectories, at each space-time we compute the 
effective mass m 2 ff- The phenomenon of long wavelength DCC amplification will occur whenever the effective mass 
squared is negative. To single out the trajectories for which maximum instabilities occur we calculate the following 
quantity, 

G = y \m e f /|6(— m 2 ^ ^rdrdxdy (13) 

This can be a measure of instability in a particular evolution. We call it amplification factor. This is an impor- 
tant parameter, as it directly relates to the size of DCC domains. In fig.l, we have shown the distribution of the 
amplification factor G for the 500 trajectories. The distribution is more or less Gaussian like. All the trajectories 
shows appreciable instability. The maximum and the minimum instabilities differing by ~ 10%. Unlike 1-dimensional 
expansion, DCC domain formation is a distinct possibility in 2-dimensional expansion scenario. It is interesting to 
note that presence of friction and noise do not change the scenario to a great extent. 

In fig. 2, we show the contour plot of 7r 2 field for the most unstable orbit. Initially, at 1 fm, the field arc randomly 
distributed. There is no domain like structure. After a few fermi of evolution, correlation starts to build up, and a 
domain like structure can be clearly seen. The domain like structure is most distinguished at 5 fm, where we see clear 
two domain formation. It may be noted that the two domains have different sign. The domain like structure persist 
at 7 fm also. The result is a clear indication of DCC like phenomena, even in presence of dissipation and noise. 

In fig. 3, we have shown the time evolution of the correlation function, defined as, 

IVrW x 

where the sum is taken over those grid points i and j such that the distance between points is r. At 1 fm, there is 
no correlation beyond the lattice spacing (.25). A long range correlation emerges at 3 fm. The correlation length 
increases as the evolution time increases and at 5 fm it is as large as 5 fm, corroborating the pictorial result of fig. 2 

To summarize, we have investigated the effect of friction on the possible DCC domain formation. In the equation 
of motion for the linear sigma model fields, we include a friction and a white noise term and solve it assuming 
boost invariance. The noise term is required to be consistent with the fluctuation-dissipation theorem. Initial field 
configuration was assumed to be Gaussian random with quenched condition (< 4> >= < <j> > = 0). At each space- 
time rn 2 ff was calculated. Phenomenon of long wavelength amplification occurs when it becomes negative. For each 
trajectories, we calculate the amplification factor, as defined in eq.Q. Amplification factor gives an indication of 
instability in the system. It was seen that for the 500 trajectories, the maximum and minimum instabilities do not 
differ much, indicating possibility of domain like structure formation in all the events. Clear domain like formation 
is seen in the event with maximum instability. For initial random field distribution, domain like structure emerges 
after a few fermi of evolution. Most distinct domain like structure is seen at 5 fm. There domain of 5 fm size is 
seen to be formed. We have also calculated the correlation function for the pion field. Similar result is also obtained 
there. Initially correlation donot extend beyond .25 fm, (the lattice spacing). After a few fm of evolution, long range 
correlation starts to build up. At 5 fm, correlation length of 5 fm is obtained. The result clearly indicate that the 
DCC domain formation is a distinct possibility, even in presence of friction and noise. 
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FIG. 1. Distribution of the amplification factor G for the 500 trajectories calculated. 
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FIG. 2. Contour plot of evolution of 7T2 field for the event with maximum instability. At 5 fm, two domain like structure can 
be seen clearly. They have opposite sign. 



6 



Correlation 



o.o 



W\ 

v \ 
\ \ 
\ \\ 

\ 


ii ii ii ii i 


1 fm 
3fm 
5fm 
7fm 









2 4 6 8 10 
r(fm) 

FIG. 3. Evolution of the correlation function. Initially at 1 fm, there is no correlation beyond .25 fm, the lattice spacing. 
Large correlation devolop at later times. 
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